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Abstract 

In this article, we prove that exact representations of dimer and plaquette valence-bond ket 
ground states for quantum Heisenberg antiferromagnets may be formed via the usual coupled 
cluster method (CCM) from independent-spin product (e.g. Neel) model states. We show that we 
are able to provide good results for both the ground-state energy and the sublattice magnetization 
for dimer and plaquette valence-bond phases within the CCM. As a first example, we investigate 
the spin-half J1-J2 model for the linear chain, and we show that we are able to reproduce exactly 
the dimerized ground (ket) state at J2/J1 = 0.5. The dimerized phase is stable over a range of 
values for J2/J1 around 0.5, and results for the ground-state energies are in good agreement with 
the results of exact diagonalizations of finite-length chains in this regime. We present evidence of 
symmetry breaking by considering the ket- and bra-state correlation coefficients as a function of 
J2I J\- A radical change is also observed in the behavior of the CCM sublattice magnetization as 
we enter the dimerized phase. We then consider the Shastry-Sutherland model and demonstrate 
that the CCM can span the correct ground states in both the Neel and the dimerized phases. 
Once again, very good results for the ground-state energies are obtained. We find CCM critical 
points of the bra-state equations that are in agreement with the known phase transition point for 
this model. The results for the sublattice magnetization remain near to the "true" value of zero 
over much of the dimerized regime, although they diverge exactly at the critical point. Finally, we 
consider a spin-half system with nearest-neighbor bonds for an underlying lattice corresponding to 
the magnetic material CaV409 (CAVO). We show that we are able to provide excellent results for 
the ground-state energy in each of the plaquette-ordered, Neel-ordered, and dimerized regimes of 
this model. The exact plaquette and dimer ground states are reproduced by the CCM ket state in 
their relevant limits. Furthermore, we estimate the range over which the Neel order is stable, and 
we find the CCM result is in reasonable agreement with the results obtained by other methods. 
Our new approach has the dual advantages that it is simple to implement and that existing CCM 
codes for independent-spin product model states may be used from the outset. Furthermore, it 
also greatly extends the range of applicability to which the CCM may be applied. We believe that 
the CCM now provides an excellent choice of method for the study of systems with valence-bond 
quantum ground states. 
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I. INTRODUCTION 



Lattice quantum spin models not only provide useful models of many physically realizable 
magnetic systems but also serve as prototypical models of strongly interacting quantum 
many-body systems. Indeed, the basic models of quantum magnets are given by lattice spin 
models that often display rich quantum phase transitions between ground states of different 
order as some control parameter is varied. Their collective behavior is extremely complex 
due to the presence of strong quantum effects. Furthermore, the underlying crystallographic 
lattices for these materials may exhibit complex symmetries. Their rich phase diagrams 
include exotic phases of novel quantum order due to the strong interplay between competing 
interactions and large quantum fluctuations. For all of these reasons they have naturally 
provided an excellent test-bed where the various methods of quantum many-body theory 
can be applied and further refined. 

Of particular interest is the formation of dimer- and plaquette-ordered singlet ground 
states (so-called valence-bond crystal (VBC) states) in quantum spin systems. Often, the 
formation of enhanced dimer or plaquette correlations is driven by frustration, which can 



increase quantum fluctuations and which may result in such 



quantum paramagnetic states 
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24| . Usually, VBC states are complicated quantum many-body states, see, e.g. 



the Heisenberg antiferromagnet on the star lattice [9|, |23|, |24|. However, for certain systems 
the VBC states are simple exact product eigenstates of the underlying Heisenberg interaction 
Hamiltonian. Examples for the appearance of such exact VBC product eigenstates are the 
spin-half J1-J2 model on the linear chain at the point J2JJ1 = 0.5 (the 



so-called Majumdar-Ghosh 
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and the Shastry- Sutherland model 
22l |. Furthermore, it is often useful to distinguish between 
VBC phases that have the same translational symmetry as the spin Hamiltonian and those 
that spontaneously break the symmetry of the underlying spin lattice. Examples of the 
former case are the Shastry-Sutherland model and the Heisenberg antiferromagnet on the 
star lattice, whereas the J1-J2 model on the linear chain is an example of spontaneous 
symmetry breaking. 

Another mechanism for the formation of non-magnetic dimer or plaquette VBC ground 



states that does not involve frustration is the competition between non-equivalent antifer- 
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romagnetic nearest-neighbor bonds. This may lead to the formation of local singlets of 



;wo (or four) coupled spins if the 
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IS of the non-equivalent bonds differ sufficiently 



361 ] . By contrast to frustration, which yields com- 



petition in quantum as well as in classical spin systems, this type of competition is present 
only in quantum systems. The symmetry of the ground state follows the symmetry of the 



Hamiltonian in such cases. An example of this is given 



the magnetic material CaV409 of Ref. 
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the nearest-neighbor models for 
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3l| in which one finds two 



non-equivalent antiferromagnetic nearest-neighbor bonds J and J' belonging to dimers (J') 
or to four-spin plaquettes (J). If these non-equivalent nearest-neighbor bonds are different 
in strength then local singlet formation may destroy Neel long-range order. As J' / J becomes 
large, local singlets are formed on dimers. Alternatively, local four-spin singlets on the pla- 
quettes are formed for strong plaquette bonds J. Another example studied in the literature 
is the so-called J- J' model on the square lattice, i.e., a model with a reg ular distribution 

n n n n n □ 

of two different nearest-neighbor bonds on the square lattice |32l . 1331 . I34I . 1351 . 1361 . I37| . In 
both cases, the formation of local singlets suppresses the magnetic Neel long-range order. 
However, the VBC ground state is now a complex many-body state; with a simple product 
VBC state appearing only in the limits J' / J ^ oo and J/ J' oo. 



In this article we focus on the application of the coupled cluster method (CCM) to 



quantum Heisenberg anti: 



'erromagnets 




one of the most powerful and most widely applied techniques of modern-day microscopic 
quantum many-body theory. The CCM allows the inclusion of multi-particle correlations 
into the ground- and excited-state wave functions in a controlled and systematic manner. 
It has been applied to a great variety of different lattice quantum spin systems with great 
success. In particular, it has been used successfully with model (or reference) states built 
by independent-spin product states for which the choice of state for the spin on each site is 
formally independent of the choice of all others. Often for these independent- spin product 
model states the use of collinear states, such as the Neel state, is possible where all spins 
are ali gned parallel or antiparallel to one axis (e.g., the z-axis) , see, Refs. |47l. |48U49U50U5U 

52I, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, y, yi . 

However, noncoUinear (e.g. spiral) model states can be favorable for certain frustrated 
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0]. Multi-spin correlations are then included 



systematically on top of the independent-spin product model states. 

The CCM for independent-spin product model states may be applied to 



ligh orders 



by using a computational implementation described in Refs. [55|, |59|, |63|, |68|, |71|, |72|. In 
particular, it may be applied to lattices of complex crystallographic symmetry. Furthermore, 
it is not constrained to systems with spin quantum number s = 1/2. 

In previous work, non-classical VBC ordering has also been considered using the CCM by 
emp loying directly valence-bond model states, i.e. two- or four-spin singlet product states 
53(1 • However, this earlier approach involves the direct use of products of localized states 
(e.g., two-spin dimers or four-spin plaquettes) in the model state. Hence, this approach 
requires that a whole new matrix-operator formalism be created for each new problem. 
Also, the Hamiltonian and CCM ket- and bra-state operators must be written in terms 
of this new matrix algebra. The CCM equations may be derived and solved once the 
commutation relationships between the operators have been established. Although formally 
straightforward, this process can be tedious and time-consuming. Furthermore, the existing 
high-order CCM formalism and codes also need to be amended extensively for each separate 
model considered. 

In this article we use a quite different way to describe VBC states. Starting directly 
from collinear independent-spin product model states, we discuss how we can form exact 
local dimer or plaquette ground states within the CCM. This approach has the advantages 
of being conceptually simple and thus also of being easy to implement. Furthermore, we 
can use directly the existing high-order CCM formalism, computer codes, and extrapolation 
schemes. After first describing the new methodology, we apply the method to specific spin 
models considered in the literature that exhibit dimer and plaquette ground states. By 
comparison with existing results, we demonstrate that the CCM provides good results for 
the ground-state properties of these systems. We then conclude the article by considering 
the implications of our results. 
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II. METHOD 



A. CCM Formalism 

The exact ket and bra ground-state energy eigenvectors, |\[^) and (\&|, of a general many- 
body system described by a Hamiltonian H, 

H\m)=Eg\^)- {^>\H = Eg{^>\ , (1) 

are parametrized within the normal (NCCM) version of the single-reference CCM as follows: 

|vl>) = e^|$) ; S = Y.SiCi, 

(§| = (<l>|^e-^ ; S = l + J2SjCY. (2) 

The normalized single model or reference state |$) is defined with respect to a suitable set of 
(mutually commuting) many-spin creation operators {Cf}. We note that the model states 
are generally related to the classical ground states of the lattice spin system. These states 
are products of single-spin eigenstates of some appropriately defined operator s| whose 
direction in a set of global axes can differ from site to site, e.g., a typical such states is: 
1$) = ■ ■ ■ ® I T) ® I i) ® I T) ® I i) ® ■ ■ ■ (Neel state). However, we remark also that planar 
model states or spiral model states in the global spin coordinate axes (|$) = ■ ■ ■ ® | — »• 
) ® I ^) ® I ® I ^) ® ■ ■ ■ and 1$) = ■ ■ ■ ® I t) ® I /) ® I ^) ® I \) ® I i) ® I / 
) (E) ■ ■ ■, respectively) may also be considered. In order to make the CCM calculations easier 
to carry out in practice, we generally rotate the local axes of the spins so that they all 
appear notationally to point in the downwards z-direction. The model state is then given 



by: 1$) 



Re 



56 



s. 



22 



34 



54 



67. 



69 



71 



74 



) ® I i) ® ■ ■ ■ . The interested reader is referred to 
77j ] for more details about spiral model states and Ref. 
801 ] for more details about planar model states. 
The operators C/ = (C7)^ with Cq = 1, have the property that ($|C/ = = 
C7|$) V / 7^ 0. They form a complete set of multi-spin creation operators with respect 
to the model state |$). Thus, the creation operators are represented simply as a product of 
spin-raising operators = sf. + is^ over the set of lattice sites {k} after rotation of the local 
frames such that all spins appear to point downwards, as described above. The creation 
operators are now given by C/ = s^s^ ■ ■ ■ s^. We note that the definitions of Eq. ([2]) imply 
the normalization {^l^) = ($|$) = 1. 
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The determination of the correlation coefficients {Sj, Sj} is achieved by requiring the 
ground-state energy expectation functional H{{Si,Si}) = = ($|S'e-^ife'^|$) to 

be stationary with respect to variations in each of the (independent) variables of the full 
set. We thereby derive the following coupled set of equations: 

6H/6S1 = ^ {(^\CYe-^He^\(!?) = 0, V/ ^ ; (3) 
6H/6Si = =^ ($|5e-^[iJ,C/]e^|$) = 0, V/ ^ . (4) 

Equation ([3]) also shows that the ground-state energy at the stationary point has the simple 
form 

Eg = {^\H\^) = {^\Se-^He^\^) 
= ($1(1 + ^5/Cf)e-^iJe^|$) 

= {^\e-^He^\<!>) + J2Si{^\CYe~^He^\^) 

^Eg = {^e-^He^\^) . (5) 

We see that the expectation value of the ground-state energy of Eq. ([5]) contains terms in 
S only and so it also contains ket-state correlation coefficients only. Generally, however, we 
need to use both the bra and ket states to find a ground-state expectation value. Indeed, 
the ground-state energy is the only special case that requires just the ket-state alone. We 
note also that this (bi-)variational formulation does not lead to an upper bound for Eg when 
the summations for S and S in Eq. ([2]) are truncated, due to the lack of exact hermiticity 
when such approximations are made. However, it is also important to realize that the 



Hellmann-Feynman theorem is preserved in all such approximations (46|. 

We note that any practical calculation requires an approximation to be made for both S 
and S. The three most common schemes are: (1) the SUBn scheme, in which all correlations 
involving only n or fewer spins are retained; (2) the SUBn-m sub-approximation, in which 
all n-spin-flip clusters spanning a range of no more than m adjacent lattice sites are retained; 
and (3) the "localized" LSUBm scheme, in which all multi-spin correlations over distinct 
locales on the lattice defined by m or fewer contiguous sites are retained. 

We use here an order parameter M = — | | ^) that is thus defined to be the 
negative of the ground-state expectation value of the operator = Yl!i=i Si- Here N 00 
is the total number of spins on the lattice and, very importantly, each operator is defined 



with respect to the local rotated spin axes on lattice site i, which are themselves defined by 
the choice of particular model state |$), as explained previously. Clearly, by definition, for 
the original model state, the order parameter M is simply — ($ | | $) = |. Quantum 
correlations (or fiuctuations) in the exact ground state will thus have the effect of decreasing 
M from this maximal value. For the sake of ease of use and clarity we shall henceforth refer 
to the order parameter M as the sublattice magnetization, although this terminology is 
strictly only appropriate for antiferromagnetic ordering on a bipartite lattice. Hence, the 
sublattice magnetization M in the local axes of the spins in which all spins in the model 
state 1$) point in the negative ^-direction is given by 

1 ^ ~ 

m=-^j:{^\su^). (6) 

This quantity is easily determined once the bra- and ket-state equations have been solved. 
B. Construction of CCM Valence-Bond Ket Ground States 

Following on from our discussion above of the use of independent-spin product states as 
CCM model states, we now present a method for creating VBC states within the CCM from 
such independent-spin product model states. As an example, we consider the Heisenberg 
model that has a Hamiltonian defined by 

H=Y.s, ■ s, , (7) 

(id) 

where the indices i and j run over all lattice sites on a given lattice. The brackets around 
indicate that all nearest-neighbor pairs are counted once and once only. For the bipartite 
lattices considered here (namely, the linear chain, the square lattice, and the "CAVO" 
lattice), we choose a model state in which nearest-neighbor spins are antiparallel along, say, 
the 2;-direction. The local frames of the "up" spins are rotated by 180° so that they point 
downwards in these local axes. This is achieved by carrying out the following transformation 
of the local axes of these spins: — > —s^; s^; and, —s^. The model state is now 

formed from a product purely of "down" spins in the rotated spin coordinates, as described 
above. 

The Heisenberg Hamiltonian is now written in terms of the new spin axes by, 

H = -T.i^ls^ + \{^i^] + 4stY) ■ (8) 

(id) ^ 
8 



We note that the total spin in the "global" ^r-direction, sf. = J^i^ s^, is a conserved quantity 
of the ground state in all of the models studied here. 

We now present the method for creating VBC states from independent-spin product 
model states. Let us consider for a moment the one-dimensional spin-half Ji - J? Heisenberg 
antiferromagnet. The relevant model state is the coUinear Neel state 0, [m, 77, 78|. However, 
as mentioned above, the ground state for J2/J1 near to 0.5 breaks the translational lattice 
symmetry and is two-fold degenerate (see also Sect. IIII Ap . In order to take into account 
this property, we have to double the unit cell, i.e. the relevant unit cell has two neighboring 
sites at points (0,0,0) and (1,0,0) and a single Bravais vector of (2,0,0)"^. Note, however, that 
such an explicit increase of the unit cells is not necessary for VBC phases that do not break 
the translational lattice symmetry. The doubling of the unit cell, now enables us to consider 
two distinct types of two-spin nearest-neighbor ket-state correlation coefficients; in this way 
allowing to break the translational lattice symmetry. We call the two nearest-neighbor 
ket-state correlation coefficients S2 and where, we define iSf to connect those nearest- 
neighbor sites between different unit cells and S2 to connect those nearest-neighbor sites 
within each unit cell. With respect to the rotated spin coordinates, we may now construct 
via Eq. ([2]) a simple dimerized product CCM ket state given by either iSg = 1 and 1S2 = 
or 5^ = and S!^ = 1 for S2 = Ei, 4a4a+i + ^iT^i^ ^heie ia runs over all 

sites with odd-numbered indices and % runs over all sites with even-numbered indices. It is 
obvious, that this choice: (i) breaks the lattice symmetry; and, (ii) represents two different 
degenerate states. The proof that the above choice for 5f , S2 leads to dimerized product 
states is as follows: 

I*) = e^2|$^ 

= (1 + .+4)(i + 44)(i + 44) ■ ■ ■ 1$) 

= {|ili2)+|TlT2)}®{|i3i4)+|T3T4)} 

® • • • (9) 

where the notation | J,j) and | ti) indicates a 'down' or 'up' spin, respectively, localized to 
site i. We note that {sf)'^\^) = V i is assumed in Eq. (Q, which holds true for spin-half 
systems such as those considered here. 



Furthermore, we note that if the local axes of spins on one sublattice are "re-rotated" 
such that 'down' spins become 'up' spins once again, i.e., so as to regain the Heisenberg 
Hamiltonian of Eq. ([7]) in the global spin coordinate system, then the dimerized product 
state becomes the usual product of the nearest-neighbor dimer singlets, namely, (| 
) I iiti)) at nearest-neighboring sites i and j. Note that if such a dimerized product state 
becomes the true ground state of a given spin problem (e.g., the spin-half linear chain J1-J2 
model at J2/Ji=0.5 - see below for more details) then for any level of LSUBm (m > 2) 
approximation either 1S2 = 1 with all other ket-state correlation coefficients (i.e., including 
those for m spins with m > 2) equal to zero or 1S2 = 1 with all other ket-state correlation 
coefficients equal to zero, is a valid ket ground-state solution to this problem. 

The situation is different for the bra state. Firstly, we note that the dimerized product 
bra state being equivalent to the corresponding dimerized ket product state of Eq. iQ in 
rotated coordinates ought to be, 

(^1 = {(ili2| + (TlT2|}®{(i3i4| + (T3T4|}®--- 

= {(ili2i3i4i5i6 ■ ■ ■ I + (TlT2i3i4i5i6 " ' ' | + (ili2T3T4i5i6 ' " ' | + 

• ■ ■ +(iii2T3T4T5T6T7T8 • ■ ■ I + ■ ■ ■ + (T1T2T3T4T5T6T7T8 ■■■!}. (10) 

We notice now that the modes of action of the spin operators (leftwards) on the bra spin 
states for s = 1/2 are, 

(ii,+ = o ; ai.- = (Ti 

(T|s- = ; (T|s+ = (i|. (11) 

Thus, the NCCM ground bra state for LSUBm with m a finite number can only ever contain 
a maximum of m "up" states in the bra state because of the linear nature of the bra-state 
operator 5* in Eq. ([2]). Hence, by contrast to the ket state, within the LSUBm approximation 
we can never construct an equivalent simple dimerized product bra state using the NCCM 
except in the exact limit where m ^ 00. However, carrying out CCM calculations in the 
limit m — s> 00 using computational methods is generally not practical. We note that this 
problem might be alleviated by using the extended coupled cluster method (ECCM). In 
this method, the bra state is written in terms of an exponential with respect to both the 
ket-state and bra-state correlation operators. It is this exponential term that allows such 
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FIG. 1: COM results at the LSUB12 level of approximation for the ground-state nearest-neighbor 
ket-state correlation coefficients of the spin-half J1-J2 antiferromagnet on the linear chain. The 
nearest-neighbor coefficients ^2 and 5^ of the symmetry breaking dimerized solution are shown 
by dashed lines. Results for the usual ('Neel-type') solution (where ^2 = are shown by the 
full line. Below the (bifurcation) COM critical point at J2I J\\c\ there is only the solution with 
^2 = S\. A termination point J2I J\\t of the CCM equations for the dimerized solution, at which 
point the real solution to the CCM equations terminates, is indicated by the boxes. 

dimer solutions, e.g., for the ket state for the NCCM above. However, we do not discuss the 
ECCM further in this present paper. 

We mention that the sublattice magnetization M (see Eq. ([H])) for the simple dimerized 
product ket state (i.e., iSg = 1 and all other coefficients equal to zero) is written in terms of 
only one bra-state correlation coefficient, i.e. M = | — IS^- 



III. RESULTS 



A. The Spin-Half Ji J2 Model on the Linear Chain 

The Hamiltonian for this spin-half model has nearest-neighbor bonds of strength J\ and 
next-nearest-neighbor bonds of strength J 2- We use a Neel model state in which nearest- 
neighbor spins on the linear chain are antiparallel. We rotate the spin coordinates of the 
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FIG. 2: CCM results at the LSUB12 level of approximation for the ground-state nearest-neighbor 
bra-state correlation coefficients of the spin-half J1-J2 antiferromagnet on the linear chain. The 
nearest-neighbor coefficients iSf and iSl of the dimerized solution are shown by dashed lines. Results 
for the usual ('Neel-type') solution (where S2 = 51) shown by the full line. Below the critical 
point at J2/Ji\ci both solutions coincide. Results for these bra-state correlation coefficients diverge 
at the critical point J2/Ji\ci- A termination point at J2/Ji\t is shown by the boxes on the right- 
hand side of the figure. 

'up' spins so that notationally they become 'down' spins in these locally defined axes. The 
relevant Hamiltonian in rotated coordinates is then given by 



H = -Ji E(^'^.' + li^TsJ + sfsj}) + J2 E + l{4s, 

{id) im) 



Si 4}) 



;i2) 



where runs over all nearest-neighbor sites on the lattice counting each pair once and 
once only and {{i, k)) runs over all next-nearest-neighbor sites on the lattice, again counting 
each pair once and once only. Henceforth we put Ji = 1 and consider J2 > 0. 

The ground-state properties of this system have been studied using methods such as 



exact diagonalizations 
approaches 



54, CCM (5]J 



53 



i?i DMRG [ai, y, y, 

5| (see Refs. 13, B for a general review). Note that previous CCM studies 



771 ] ■ and field-theoretical 



of the model considering only independent-spin product model states that conserve the 



lattice symmetry are reported in Refs. 



51 



78|. At J2/-^i = we have the unfrustrated 



Heisenberg antiferromagnet, where the exact solution is provided by the Bethe Ansatz. The 
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FIG. 3: CCM results for the ground-state energy of the spin-half J1-J2 antiferromagnet with Ji = 1 
on the linear chain at the LSUB12 level of approximation. The dimerized and usual ('Neel-type') 
solutions are shown in this figure. Results of exact diagonalizations for = 28 and = 32 number 
of sites is also shown. The CCM termination point for the dimerized solution is shown by the box. 

ground state is gapless and the spin-spin correlation function (sj ■ s^) decays slowly to zero 
according to a power-law, i.e. no true Neel-like long-range order is observed. In the region 
J2/J1 > the nearest-neighbor (Ji) and next-nearest-neighbor interactions (J2) compete, 
thus leading to frustration. At J2I J\ = 0.2411(1) the model exhibits a transition to a two- 
fold degenerate gapped dimerized phase with an exponential decay of the correlation function 

nnnri 

(Sj ■ Sj) [2|, y, |5|, l6[. This state breaks the translational lattice symmetry. At the Majumdar- 
Ghosh point J2/ Ji = 0.5 there are two degenerate simple exact dimer-singlet product ground 
states corresponding to the dimerized product state of Eq. iQ for the Hamiltonian of Eq. 
( IT2|) [l|. (We recall that rotated spin coordinates are used in Eq. ([9]).) 

We now consider how this model can be treated in the dimerized phase by the CCM via, 
as presented above, the identification of a special dimerized solution of the CCM equations 
for a Neel model state. Also as discussed above, we use a doubled unit cell including two 
neighboring sites for a spin-half system on the linear chain at points (0,0,0) and (1,0,0) 
and a single Bravais vector (2,0,0)"^ to take into account the symmetry breaking. There 
are thus two distinct types of two-spin nearest-neighbor ket-state correlation coefficients 
and again these are denoted as 5f and S2. The exact ground state at J2/J1 = 0.5 of Eq. 
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FIG. 4: CCM results for the sublattice magnetization M of the spin-half J1-J2 antiferromagnet 
on the linear chain. Below the critical point at J2/J1IC1 the results for the usual ('Neel-type') and 
the dimerized solution coincide. At J2/J1IC1 the sublattice magnetization of the dimerized solution 
exhibits a jump, whereas M for the usual ('Neel-type') solution is continuous. Above J2/Ji\ci the 
sublattice magnetization of the dimerized solution (black lines) is smaller than that of the usual 
('Neel-type') solution (grey lines), except for the LSUB8 in a small region above J2/J1 = 0.5. 

([9]) is obtained by setting 1S2 = 1 and all other coefficients equal to zero. Starting from 
J2/ Ji = 0.5 we are able to track this exact solution at J2/ Ji = 0.5 within a certain LSUBm 
approximation for other values of J2/J1 and the results for the nearest-neighbor ket-state 
correlation coefficients in LSUB12 approximation are presented in Fig. [TJ Clearly, we see 
that the exact dimerized product-state solution for the ket ground state is obtained within 
LSUB12 level of approximation (and, indeed, at all LSUBm approximation with m > 2) 
at J2/J1 = 0.5, i.e. ^2 = 1 and all other coefficients equal to zero. Moving away from 
J2/J1 = 0.5 we still find a CCM ground state that breaks the lattice symmetry. However, 
this dimerized state deviates from the simple product state, i.e. ^2 7^ 1 and other non-zero 
coefficients Sj occur. Furthermore, the solution (i.e. S2 = S2) having full translational 
symmetry is the only solution below a critical point J2/Ji\cii< 0.5). Henceforth, we shall 
refer to this solution below J2/Ji\ci as the "usual ('Neel-type') solution" because previous 



CCM calculations 
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531] for the J1-J2 model have considered the non-symmetry breaking 
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case only. For larger values of J2/ Ji a CCM termination point is observed at J2/ Ji\t (> 0.5), 
shown by the boxes in Fig. [H At this point, the real solution of the CCM dimerized solution 
is terminated. These CCM results indicate that a dimerized phase exists over a finite range of 
J2/ Ji, which is in agreement with known results, see e.g. Refs. 0,5,6]. Qualitatively similar 
results are observed at other levels of LSUBm approximation for the ket-state correlation 
coefficients as a function of J2I J\- The results for J2/ J\\cx ^-iid J2I J\\t are shown in Tabled 
It is obvious that the CCM critical point J2I J\\cx becomes smaller (i.e., becomes closer to the 
true critical point J2/ J\ = 0.2411(1) j^, ^) with higher orders m of LSUBm approximation. 
However, the critical point J2/ Ji\ci is still sigiiificantly too high even at the LSUB12 level of 
approximation. However, as shown in Ref. 5|] the dimerization gap for J2/J1 < 0.4 is very 
small and it is therefore not very surprising that we do not detect the dimerized phase below 
J2/J1 < 0.4 using LSUBm approximations with m < 12. On the side of J2/J1 > 0.5 the 
existence of termination points can be related to the appearance of incommensurate spiral 
spin correlations at J2/ Ji > 0.6 

mm 

that are not taken into account in the model 

state used here. 

The nearest-neighbor bra-state correlation coefficient at the LSUBm level of approxima- 
tion at J2/J1 = 0.5 has = 1/4 with m > 4. This is shown in Fig. [2] for the LSUB12 
level of approximation. We find that the bra-state solution for the nearest-neighbor corre- 
lation coefficients is close to 1/4 over the range J2/Ji\ci < J2/J1 < J2/Ji\t- However, we 
find that the nearest-neighbor correlation coefficient diverges as J2/J1 J2/ Ji\ci and this 
is also shown in Fig. [2l Again, the usual ('Neel-type') solution (iSf = S\) is obtained for 
J2/ Ji < J2/ Ji\ci- The upper CCM termination point at J2/ Ji\t is also shown in Fig. [2] by 
the boxes on the right-hand side of the figure. Once more, qualitatively similar results are 
observed at other levels of LSUBm approximation for the bra-state correlation coefficients 
as a function of J2I J\- 

We now consider the ground-state energy of this system in the dimerized regime. Our 
results for the new dimer solution and the usual ('Neel-type') solution are shown in Fig. [3l 
Firstly, we note that the exact ground-state energy of Eg/N = — 0.375Ji is obtained at 
the point J2/J1 = 0.5, as expected. We note again that our solution is an exact ground 
eigenstate at this point. We see that ground-state energy of the usual ('Neel-type') CCM 
solution in which ^2 = S2 at the LSUB12 level of approximation actually lies below this 
exact solution. This indicates (i) that the usual ('Neel-type') CCM solution is a relatively 
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TABLE I: CCM results for the positions of the range of the dimerized phase. 



LSUBm ,/2/JiU J2/J1 



2 


0.4761 




4 


0.4745 


0.5576 


6 


0.4637 




8 


0.4568 


0.7410 


10 


0.4498 


0.6404 


12 


0.4429 


0.5956 



poor choice at this point; and, (ii) that the CCM ground-state energy does not fulfill the 



variational principle 



6l| . Furthermore, we see that CCM dimer solution compares extremely 



well to results of exact diagonalizations for = 28 and = 32 sites in the dimerized regime 
shown in Fig. [H It certainly provides far better results than those of the usual ('Neel-type') 
CCM solution beyond the critical point at J2/ Ji\ci- 

The results for the sublattice magnetization M of this model are presented graphically in 
Fig. m Since the one-dimensional J1-J2 model does not possess Neel long-range order for 
any value of Jj, J2 > the true value is M = 0. As is known from previous CCM calculations 



51 



the sublattice magnetization is nonzero (but small) using the usual Neel 
model state. However, the correct result M = can be obtained 77|, l78[ by extrapolating 
the 'raw' LSUBm data to m ^ 00. Indeed it is obvious that the CCM-LSUBm values for M 
are non- negligible for Neel model state in the region J2/J1 < J2/ Ji\ci- It is also obvious that 
M decreases with the level of approximation m approaching the true value M = 0, and, that 
increasing the strength J2 of the frustration weakens magnetic order. More interestingly, we 
find that the sublattice magnetization behaves disco ntinuously at J2/ Ji\ci by tracking the 
lattice symmetry-breaking dimerized solution, and then remains near to zero at LSUBIO and 
LSUB12 levels of approximation across the entire range J2/ Ji\ci < J2/J1 < J2/ Ji\t- (At the 
lower LSUB8 level of approximation the results for the sublattice magnetization differ from 
zero by a small amount in a small region above J2/ Ji > 0.5.) On the other hand, by tracking 
the usual ('Neel-type') solution M changes continuously with J2 and it is larger than for 
the dimerized solution for J2/Ji\ci < J2/J1 < 0.5. This behavior of M is another indication 
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FIG. 5: The nearest-neighbor bonds (sohd hnes) of strength Ji and the next-nearest-neighbor 
diagonal bonds (dashed hnes) of strength J2 for the Shastry-Sutherland model. The geometric 
unit cell is shown by the square with the grey dotted lines. 

that the dimerized CCM solution describes the true physics of the model much better than 
the usual Neel solution. We note finally that the CCM sublattice magnetization is exactly 
zero at the Majumdar-Ghosh point J2/ Ji = 0.5 at all levels of LSUBm approximation using 
the dimerized product state. 

We remark again that our results for the lower phase transition point J2/ Ji\ci overestimate 
the position by a factor of two. We re-iterate that we believe that we over-estimate J2/ JAci 
because the energy gap only becomes large for values of J2/ Ji of approximately 0.4 |5|. 
Finally, we note that results given here present the possibility that the CCM might be applied 
to detect spontaneous symmetry breaking for systems of two or three spatial dimensions, 
i.e., where other approximate methods become less accurate or may not even be applicable 
(e.g., such as the DMRG method). 



B. The Shastry-Sutherland Antiferromagnet 

Another model that demonstrates dimer order is the Shastry-Sutherland antiferromag- 



net 



The Shastry-Sutherland antiferromagnet is a spin-half Heisenberg model on an 
underlying square lattice with antiferromagnetic nearest-neighbor bonds Ji and with one 
antiferromagnetic next-nearest-neighbor diagonal bond J2 in every second square (plaque- 
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tte), as shown in Fig. [51 We note that for bonds of equal strength, i.e., Ji = J2, the 
Shastry-Sutherland model is equivalent to a Heisenberg model on one of the eleven uni- 
form Archimedean lattices 19|. Interest in this model has been renewed by the discovery of 

nn 

the magnetic material SrCu(B03)2 [9|, 110| that can be understood in terms of the Shastry- 
Sutherland model. The ground state of this model in the limit of small frustration J2I J\ <^ 1 
and large frustration J2/J1 ^ 1 is well understood. However, the ground-state phase at in- 
termediate values of J2/J1 ~ 1 is still a matter of discussion. 

Just as in the case of the one-dimensional J1-J2 model this model also has a simple exact 
dimer-singlet product ground state in a certain parameter region. However, by contrast with 
the one-dimensional J1-J2 model the dimer-singlet product ground state of the Shastry- 
Sutherland model is built up of a product of dimer singlets located on the next-nearest J2 
bonds and does not break the translational symmetry. 

This model has been treated previous 

diagonalization of small lattices ay , series expansions , Q, QJieL 

malization group [l^, a gauge-theoretical approach 17], and the CCM 22 . \l\ . A recent 



y by Schwinger boson mean-field theory [10| , exact 



the renor- 



review can be found in Ref. [181]. We know from these studies that the physics of the 



or small J2 < Ji, i.e., we have 
that a simple dimer-singlet 



quantum model is similar to that of its classical counterpart 
semi-classical Neel long-range order. Furthermore, we know 
product state given by dimer singlets on the diagonal bonds indicated by the dashed lines 
in Fig. [5] is the quantum ground state for large J2. The energy per site of this dimer-singlet 
product state state is -Edimcr/A^ = — 3J2/8. It becomes the ground state for J2 > J| where 
J2 ~ (1.465 ± 0.025) Ji. Note that, by contrast to the one- dimensional J1-J2 model, the 
transition to the dimerized phase in the Shastry-Sutherland model is most likely of first 
order. 

The application of the CCM to this model has been discussed at length in Ref. 22] . The 



interested reader is referred to this reference for more details about both the model and the 
details of the applying the CCM to it. We note that the CCM solution for the Neel model 
state with nearest-neighbor spin antiparallel was identified and this worked well in the region 
of J2/ Ji < 1-6. This case has a Hamiltonian similar to that of Eq. ([8]), the nearest-neighbor 
bonds Ji (solid lines) and next-nearest-neighbor bonds J2 (dashed lines) run over those sites 
on the square lattice, as opposed to the linear chain for the case presented above. 

However, we will show that a simple dimerized product ket state solution to the CCM 
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equations also exists. Note firstly that we may define a collinear independent-spin product 
model state. This is a model state in which next-nearest-neighbor spins are antiparallel. 
This state is relevant for large antiferromagnetic J2, since antiparallel next-nearest-neighbor 



spins satisfy the J2 bonds. Thus, the "u 
columns (or rows), see e.g. Ref. [tq], 80, 



d" and "down" spins form alternate neighboring 
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83| . We choose the former case and call 



the corresponding model state the 'columnar model state'. We rotate the "up" spins into 
(nominally) "down" as illustrated in Sect. Ill Al although we must now also take the columnar 
form of the model state into account. The relevant Hamiltonian in the appropriate local 
axes described above is given by 

+-^1 i^i s'i + 7;{s7 + si sj }) 

(*!/ tOy) 

-J2Y.(^>l + \{^^^k+^i^k]) , (13) 

({i,k)) 

where the sum on {ix,jx) runs over all nearest-neighbor pairs of the lattice sites in the row or 
x-direction and the sum on {iy,jy) runs over all nearest-neighbor pairs of the lattice sites in 
the column or y-direction. Furthermore, the sum on {{i, k)) runs over distinct next-nearest- 
neighbor pairs of sites connected by the broken lines in Fig. |5l We count each bond once 
and once only (for both the nearest-neighbor and next-nearest-neighbor bonds). The unit 
cell for this model contains four sites, and it is shown also in Fig. [51 Once again, in what 
follows we set Ji = 1 and treat J2 > as the parameter of interest in the model. 

At J2 = 0, we have the usual square- lattice antiferro magnet, and the system is Neel- 
ordered at this point and in a finite range for J2 > 0. Furthermore, we note that the ground 
state becomes the simple dimer-singlet product state for J2 > J| where J| ^ (1.465 ± 
0.025) Ji. We are able to define CCM correlation coefficients for spin dimers located on the 
diagonal bonds of the dashed lines of Fig. [5] with respect to the columnar model state. 
By setting these ket-state correlation coefficients to unity and all other CCM multispin 
correlation coefficients to zero, we are able to form the relevant dimer-singlet product state 
that is the exact ground state in this regime. We find that this CCM dimer solution is a 
stable solution for the CCM equations for all values of J2/J1 > 0, and this is because the 
dimer-singlet product state is a true eigenstate for any values of Ji and J2. However, the 
energy of this state is low enough for it to become the ground-state energy only for large J2. 
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FIG. 6: The CCM ground-state energy for the spin-half Shastry-Sutherland model (with Ji = 1) 
using the Neel model state (LSUB4, LSUB6, LSUB8) and the columnar model state. Note that in 
the latter case the exact ground-state energy is obtained for any LSUB?7t, with m >2. 

The CCM dimer solution yields also the correct exact energy for the dimer- singlet product 
state, namely, Eg/N = — 0.375J2- The CCM ground-state energies are now shown for both 
the dimer and usual Neel results in Fig. [61 We see that there is a crossing of the Neel and 
dimer energies at J2/J1 ~ 1.48, as reported in Ref. 7o| . 

Another interesting point is that the bra-state correlation coefficients do not remain 
constant with respect to varying J2. Indeed, we find many of the bra-state correlation 
coefficients diverge at a CCM critical point as may be seen in Fig. [7] for the next-nearest- 
neighbor bra-state correlation coefficients on the diagonals (i.e., those corresponding to the 
dashed lines in Fig. [5]). This is a critical point for the CCM bra-state equations only. The ket- 
state equations clearly do not contain a similar critical point. Hence, this is a critical point 
that is "driven" by the bra state alone. The critical points occur at values for J2/ Ji of 1.059, 
1.243, and 1.397 for the SUB2-4, LSUB4, and LSUB6 levels of approximation, respectively. 
This is in agreement with the position of the phase transition point at ~ (1.465±0.025) Ji. 
We note that the CCM ket state is an exact ground eigenstate in this regime, whereas the 
bra state is not. As mentioned before, this is because the ket states and bra states are not 
explicitly constrained to be Hermitian conjugates of each other in the CCM parametrizations. 
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FIG. 7: COM results using the columnar model state for the ground-state bra-state correlation 
coefficients for those sites connecting the dimers in the spin-half Shastry-Sutherland model. 

The results for the sublattice magnetization M (with respect to the columnar model state) 
are shown in Fig. [HI We see that the values for the sublattice magnetization are negative for 
all values of J2/J1. Note that the true values for the sublattice magnetization are actually 
zero in this regime. Negative values of M might be an indication of missing magnetic long- 
range order. However, we note that although the COM results are negative, they remain 
close to zero (e.g., \M\ < 10~^ at the LSUB6 level of approximation) for J2/J1 > 2. This 
indicates that we obtain generally good results for the sublattice magnetization. However, 
we see from Fig. [8] that our results clearly become worse for J2I J\ < 2. For example, we 
note that the sublattice magnetization diverges as we approach the CCM critical point (of 
the bra-state equations only) at J2/ Ji- We note that the bra and ket states are not explicitly 
constrained to be Hermitian conjugates. Hence, the bra state does not have to be an exact 
eigenstate of the Hamiltonian even though the ket state is for the Shastry-Sutherland model. 
We note therefore that this is enough to allow the bra-state equations to become critical even 
though the ket-state equations do not; hence, the nearest- neighbour bra-state correlation 
coefficient (^2 ) diverges. The sublattice magnetization (M = | — 2S2) therefore diverges 
also. Again, this is a reflection of the critical point that is "driven" by the bra state alone. 
The Neel state with nearest-neighbor spins antiparallel is the appropriate CCM model state 
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FIG. 8: The sublattice magnetization for the spin-half Shastry-Sutherland model for the columnar 
model state. 



22] below the critical point J2 < J[ 



C. The J— J' Heisenberg antiferromagnet on the CAVO lattice 



In this section, we consider an antiferromagnetic Heisenberg model in which the basic 
geometric unit cell contains four neighboring lattice sites on the underlying crystallographic 
lattice of the magnetic material CaV409 (CAVO), shown in Fig. [91 There are two non- 
equivalent antiferromagnetic nearest-neighbor bonds J and J' belonging to dimers (J') and 
to four-spin plaquettes (J) respectively. The ground state of the quantum model depends on 
the ratio J' / J of the competing bonds. Using a unit cell as defined in Fig. [9], the plaquette 
bonds J are inside the four-site unit cell and the dimer bonds J' connect sites in different 
unit cells. We note that this model is not frustrated but the two non-equivalent nearest- 
neighbor bonds lead to a competition in the quantum system. Henceforth, we choose an 
energy scale such that J = 1. 



We note that several techniques [25 



26 



27 



28 



29 



30 



3l| suggest that the Neel-ordered 



ground state at J' / J = 1 persists over a finite range of values of J' /J around this point, 
J^yj < J'/ J < J'c^/J- The best estimates for J'^^/ J and J'c^/ J are probably provided by 
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FIG. 9: The CAVO lattice. The nearest-neighbor bonds that connect two sites on a four-site 
plaquette are shown by the sohd hnes and have a bond strength given by J (=1). The nearest- 
neighbor bonds that connect two sites on different plaquettes (dimer bonds) are shown by the 
dotted hnes and have a bond strength given by J'. The unit cell of the lattice is shown by the 
square with the grey dashed lines. 
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FIG. 10: CCM results at the LSUB8 level of approximation for the ground-state nearest-neighbor 
two-spin A and B ket-state coefficients (^Sf and 51) and the plaquette four-spin ket-state correlation 
(tS^) coefficients of the spin-half antiferromagnet on the CAVO lattice. 
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FIG. 11: COM results at the LSUB8 level of approximation for the ground-state nearest-neighbor 
two-spin A and B bra-state coefficients ^^'^ •^2) ^^'^ ttis plaquette four-spin bra-state (Sf) 
bra-state correlation coefficients of the spin-half antiferromagnet on the CAVO lattice. 



quantum Monte Carlo calculations 



26 



271] that suggest that the range over which Neel order 



is stable is given by J'^JJ ^ 0.939 and J'^JJ ^ 1.68 ± 0.14. For J'/ J < J'^J J the model 
exhibits a quadrumerized plaquette VBC phase with enhanced spin correlations on the four- 
spin plaquettes, and for J' / J > J'^J J it has a dimerized VBC phase with enhanced spin 
correlations on the dimers. Neither of the VBC phases breaks the translational symmetry 
of the lattice. Furthermore, we mention that by contrast with the previously considered 
models (see Sees. IIII Al and IIII Bp this model has no simple exact product ground state for 
any value of J'/ J. The interested reader can find more information on the ground-state 



phases in Refs. [25, 



26|, 



27 



28 



29 



30 



3l|. 



The four-site plaquettes in the unit cell become decoupled only in the limit J'/ J = 0. 
The ground state is a product of such four-site plaquette singlets in this limit. To model 
such a state using the CCM we start again from the Neel model state; namely, a state in 
which the spins on nearest-neighbor sites are antiparallel. 

To create an exact plaquette-singlet product VBC ground state at J'/ J = using the 
CCM we have to adjust the nearest-neighbor correlation coefficients S2 and S\ and a single 
four-body plaquette correlation coefficient ^4 containing all four sites properly. (Note that 
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FIG. 12: CCM results for the ground-state energy of the J-J' Heisenberg antiferromagnet on the 
CAVO lattice (with J = 1). 

1S2 represents those ket-state coefficients for the nearest-neighbor two-body cluster connect- 
ing sites on a plaquette indicated by the sohd fines in Fig. [9|, whereas S\ represents those 
ket-state coefficients for the nearest-neighbor two-body cluster connecting sites on a dimer 
indicated by the dotted lines in the same figure. The coefficient 1S4 represents those ket-state 
coefficients for the four-body cluster corresponding to a plaquette indicated by the solid lines 
in Fig. [9l) Indeed, it is easy to show that setting the ket-state correlation coefficients iSf and 
^4 to a value of 0.5 and all other ket-state correlation coefficients (including S\) to zero the 
plaquette-singlet product VBC state is obtained exactly, see Fig. [TOl We are able to track 
this plaquette solution as J' / J is increased away from the point J' / J = where it is exact. 
Furthermore, we are also able to reproduce exactly the dimer-singlet product ground state 
in the limit J'/ J 00. In this limit, the nearest-neighbor ket-state correlation coefficient 5| 
on the dimer bonds (dotted lines in Fig. [9]) has a value of one and all other coefficients (e.g., 
iSf and 1S4) are zero, see Fig. [TOl This solution is produced automatically when we track the 
solution (outlined above) from J'/ J = 0, and so our CCM Ansatz produces accurate results 
in all phases of this model. The corresponding bra-state correlation coefficients iSg , S2 and 
1S4 behave smoothly in the entire range of J'/ J, see Fig. [TTl 
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FIG. 13: COM results for the ground-state sublattice magnetization M of the J-J' Heisenberg 
antiferromagnet on the CAVO lattice. 

Results for the ground-state energy are shown in Fig. [T2i We see that the exact ground- 
state energy at J' / J = is reproduced for LSUBm levels of approximation with m > 4, as 
expected. At J' / J = 1, we reproduce a previous result 73| using the Neel model state for this 
CAVO lattice. At all values of J' / J, the LSUBm results are seen to converge rapidly with 
increasing levels of LSUBm approximation. We note that the system should again decouple 
into dimers as J'/ J oo and the correct ground-state energy (—0.375 J') is reproduced in 
this limit. The CCM provides excellent results for the ground-state energy for all values of 
J' /J. 

The results for the sublattice magnetization are shown in Fig. [131 We extrapolate the 
raw LSUBm data to the limit m — oo in order to determine the quantum phase transition 
points where the magnetic Neel long-range order vanishes. An appropriate extrapolation 
rule for the mag netic order parameter for systems showing a ground-state order- disorder 



transition is 



78 



79 



80 



82 



M{n) = bo + + 62(1/^^)^^^, where we use 

LSUBm results with m = {4,6,8,10}. The results for these quantum critical points are 
shown in Table [TTl Again, we note that these results indicate that the Neel order persists 
over a finite range around J'/ J = 1 for the of the J-J' Heisenberg antiferromagnet on 
the CAVO lattice. Finally, these results show that the CCM can handle plaquette VBC 
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TABLE II: Results for the quantum critical points of the J -J' Heisenberg antiferromagnet on the 
CAVO lattice. 



Method 


J' IJ 


J' IJ 


CCM Extrapolation [m = {4, 6, 8, 10}) 


0.82 


1.82 


QMC [26, 27] 


0.939 


1.68±0.14 


Cumulant Series Expansions [291 


0.9 


1.7 


Non-Linear Spin- Wave Theory [311 


0.90 


1.6 


Schwinger-Boson Mean-Field Theory [30] 


0.6 


2.4 



ordering as easily as dimer VBC ordering or the usual Neel ordering. Furthermore, we have 
demonstrated that quantum critical points can be determined by using the CCM to high 
orders of LSUBm approximation. 

IV. CONCLUSIONS 

We have shown in this article that we can easily form dimer and plaquette VBC ground 
states using the CCM with independent-spin product model states. We have investigated 
a number of cases in which the ground state was a product of localized dimer or plaquette 
singlets. Firstly, we considered the spin-half J1-J2 model for the linear chain. We showed 
that we are able to reproduce exactly the dimerized ground state at J2/ Ji = 0.5. Inter- 
estingly, a symmetry-breaking dimerized CCM solution is observed for J2/J1 < 0.5, which 
only becomes equal to the usual ('Neel-type') solution that conserves the lattice symmetry 
at a CCM critical point J2/Ji\ci- Results for the bra state correlation coefficients diverged 
at this point also. We took this to indicate the onset of the dimerized ground-state phase 
that breaks the translational lattice symmetry. We found that the dimerized phase extends 
over a finite range of values of J2/J1 both above and below 0.5. Results for the ground-state 
energy for the dimerized CCM solution were found to agree extremely well with the results 
of exact diagonalizations for = 28 and = 32 chains in the dimerized regime. The 
change from the usual ('Neel-type') solution to the dimerized solution was also observed in 
the behavior of the sublattice magnetization. 
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We then considered the Shastry- Sutherland model and we demonstrated that the CCM 
can span the correct ground states in both the Neel and the dimerized phases. We found a 
CCM critical point for the dimerized solution that is "driven" by the CCM bra state alone. 
Furthermore, the position of this critical point is in good agreement with the known value for 
the phase transition point of this model. Results for the sublattice magnetization (that ought 
to be zero in the dimerized phase) were found to be generally good in the dimerized regime. 
For example, results for the sublattice magnetization at the LSUB6 level of approximation 
were found to be \M\ < 10^^ for J2/J1 > 2. However, the sublattice magnetization was also 
found to be negative and to diverge at the critical point of the bra-state equations. 

Finally, we considered a spin-half Heisenberg antiferromagnet with nearest-neighbor 
bonds with respect to an underlying lattice that corresponds to that of the magnetic material 
CaV409. The four nearest-neighbor bonds that connect sites on a four-site plaquette have 
a bond strength given by J and the nearest-neighbor bonds that connect sites belonging to 
different plaquettes (dimer bonds) have a bond strength given by J'. The exact plaquette- 
singlet ground state at J'/ J = and the exact dimer-singlet ground state as J'/ J ^ 00 
were both reproduced exactly using the CCM with the same choice of Neel model state. 
We found that the CCM can provide precise results for the ground-state energy over all 
intervening values of the parameter J'/ J. Results for the sublattice magnetization were 
presented, and these results indicated that the Neel-ordered regime persists over a finite 
range of values of J'/ J around the point J' / J = 1. For large (and small) values of J' / J, the 
Neel long-range order is destroyed by local singlet formation on dimers (and on plaquettes). 
Extrapolations of LSUBm data suggest that the Neel-order regime extends over the range 
0.82 < J' /J < 1.82. These results were found to be in fairly reasonable agreement with 



quantum Monte Carlo results for this model 



3, 



271]. However, a discussion of the accuracy 



of phase transition points estimated using the CCM is beyond the scope of this article and 
so will form the contents of another article, although we note here that we believe also that 
higher orders of LSUBm approximation would provide closer agreement. 

As noted above, the CCM is one of the most powerful and most widely applied techniques 
of quantum many-body theory. One of the reasons for this success is based on the fact that 
the CCM allows the inclusion of multi-particle correlations into the ground- and excited- 
state wave functions in a controlled and systematic manner. The range of applicability of the 
CCM to lattice quantum spin systems has been greatly extended previously by the creation 
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of efficient and powerful high-order computer codes for independent-spin product (e.g. Neel) 
model states. These codes are simple to use and they are generally accurate in practical 
applications. Furthermore, they are extremely flexible in terms of defining and solving new 
spin problems. 

Previously however, non-classical orderings (such as local singlet formation) have also 
been considered using the CCM by employing non-Neel model states. This typically in- 
volved the direct use of products of, e.g., local dimer singlets, in the model state. However, 
this approach required a whole new matrix-operator formalism to be created for each new 
problem 53|]. This is usually tedious and time-consuming, although it is normally straight- 
forward mathematically. More importantly, however, the existing high-order CCM codes 
would then need to be amended extensively also in order to implement the new matrix al- 
gebra for each new problem. Here we have presented a much simpler and more universal 
approach that combines exact solutions for dimer or plaquette VBC product grou nd states 



with the computational implementation described in Refs. [55|, l59|, l63|, l68|, O, l72] based on 
independent-spin product model states. 

One seeming shortcoming of this new approach was found to be that the ket state can 
be an exact representation of the true ground state, whereas the bra state might not be at 
the same level of LSUBm approximation. This is due to the simple fact that the NCCM 
parametrizations of the ket and bra wave functions are not manifestly Hermitian conjugates 
of each other. This meant, for example, that the exact ground-state energy of the Shastry- 
Sutherland model in the dimerized phase was reproduced, since it is calculated from the 
ket-state correlation coefficients alone, whereas the exact sublattice magnetization (known 
to be zero in this regime) was not, since its calculation requires the use of the bra-state 
correlation coefficients as well as the ket-state coefficients. We speculate that this problem 
might be overcome by employing the extended coupled cluster method (ECCM) that contains 
an exponentiated form of the bra-state correlation operator in the (ground) bra state. We 
believe that it is possible to construct a similar exact bra-state solution for the ECCM as we 
have constructed here for the ket state using the NCCM based on independent-spin product 
model states. 

As noted above, an alternative approach is to utilize model states that are formed directly 
from products of local dimer or plaquette states (even for the NCCM). In this case, the ket 
and bra state can become Hermitian conjugates (trivially) when the model state is the exact 
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ground state. Indeed, in this case, all of the CCM ket- and bra-state correlation coefficients 
become zero. However, even here, we note that the lack of manifest Hermiticity is a general 
feature of the CCM, i.e., one that can exist even for results generated by calculations based 
on such valence-bond model states. Indeed, in those cases for which the model state is not 
(trivially) the ground state, the bra and ket states are again not constrained manifestly to 
be Hermitian any given level of LSUBm approximation. 

In conclusion, this new approach for dimer- and plaquette-ordered ground states is flex- 
ible, simple to implement, and very powerful. Lattices of arbitrary complexity can also 
be treated using this new method. Furthermore, this approach is simple because we are 
using generally independent-spin product model states derived from classical ground states. 
Indeed, this is far simpler than the alternative of creating a whole new matrix formalism 
for each new model state formed from products of localized states. Finally, this approach 
is powerful because the high-order codes based on independent-spin product model states, 
which have been employed previously with great success, may be used directly in order to 
simulate the properties of these non-Neel states. The useful LSUBm and SUBn approxima- 
tion schemes devised for the Neel model states may be used directly also. The results of the 
LSUBm scheme may be extrapolated easily to the limit m oo using existing 'heuristic' 
extrapolation schemes. The results presented here offer a great enhancement to the range of 
applicability of the CCM for lattice quantum spin systems that demonstrate 'novel states' 
of quantum order. 
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